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The requirements of conformal invariance for the two point function of the energy 
momentum tensor in the neighbourhood of a plane boundary are investigated, restricting 
the conformal group to those transformations leaving the boundary invariant. It is shown 
that the general solution may contain an arbitrary function of a single conformally invariant 
variable v, except in dimension 2. The functional dependence on v is determined for 
free scalar and fermion fields in arbitrary dimension d and also to leading order in the 
e expansion about d = 4 for the non Gaussian fixed point in ^ 4 theory. The two point 
correlation function of the energy momentum tensor and a scalar field is also shown to have 
a unique expression in terms of v and the overall coefficient is determined by the operator 
product expansion. The energy momentum tensor on a general curved manifold is further 
discussed by considering variations of the metric. In the presence of a boundary this 
procedure naturally defines extra boundary operators. By considering diffeomorphisms 
these are related to components of the energy momentum tensor on the boundary. The 
implications of Weyl invariance in this framework are also derived. 



1 Introduction 

Besides the usual bulk critical phenomena in a statistical mechanical system there 
are theoretically interesting and experimentally measurable effects associated with the 
presence of a surface, or the consequences of finite size, which are inevitably present in any 
realisable physical system. There are possible phase transitions corresponding to surface 
ordering such that, at or close to a critical point where the correlation length is large 
and microscopic details are unimportant, there are additional critical exponents describing 
the behaviour of correlation functions of operators on or in the neighbourhood of the 
boundary surface. These may also be discussed within the framework of continuum scale 
invariant quantum field theories just as for conventional considerations of bulk critical 
phenomena [1] . Furthermore, as Cardy has shown [2] , it is similarly natural to impose also 
the consequences of conformal invariance which is present at the critical point assuming 
the trace of the energy momentum tensor vanishes. This restricts the functional form of 
correlation functions in the neighbourhood of the boundary although the constraints are 
less than in the bulk since, assuming a plane boundary for <i-dimensional Euclidean space 
with d > 2, the symmetry group which leaves the boundary invariant is reduced from the 
full conformal group 0(d + 1, 1) to 0(d, 1). Hence in general even two point functions of 
primary operators are not determined uniquely up to a constant factor but may contain 
an arbitrary function. 

Nevertheless Cardy also subsequently obtained significant relations [3] between criti- 
cal exponents and the coefficients of the universal terms in the Casimir energy for simple 
geometries. In part his arguments depended on the assumption that the functional depen- 
dence of the two point function for the energy momentum tensor was essentially unique 
for conformal theories in the neighbourhood of a boundary for general d. This derivation 
has recently been criticised [4] since the required assumption was incompatible with the 
results of explicit calculations to 0(e) at the non Gaussian fixed point present in scalar (p 4 
theory when d = 4 — e. 

In this paper we re-analyse the implications of conformal invariance for two point 
functions involving the energy momentum tensor in the presence of a plane boundary. 
Although the conditions we derive are necessary rather than sufficient they allow for an 
arbitrary function to be present in the two point function, except when d = 2, and this 
freedom appears to be required to accommodate the results of calculations in specific 
conformal field theories for d > 2 (even for free scalar fields the two point function depends 
on the particular boundary conditions, compatible with conformal invariance, obeyed by 
the scalar field). 

In the next section we therefore derive the consequences of conformal invariance for 
the two point function of the energy momentum tensor near a boundary by obtaining 
differential equations for the dependence on a variable v, formed from the spatial arguments 
x, x' and which is an invariant under conformal transformations leaving the boundary 
invariant. For either x or x' on the boundary v = 1 whereas for both x, x' at large 
distances from the boundary compared with their separation v — > 0. In this limit the two 
point function may be related to that for the bulk critical system neglecting any boundary 
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effects. For d = 2 the equations have a unique functional solution but not if d > 2. We 
also consider the two point function of the energy momentum tensor and a scalar field of 
arbitrary dimension which may be non zero in the presence of a boundary and for which 
there is a unique functional solution for the dependence on v for any d. In section 3 
we determine the functional dependence on v for free scalar fields, with either Dirichlet 
or Neumann boundary conditions which maintain conformal invariance, corresponding to 
the simple Gaussian fixed point in scalar field theories for any d. We also calculate the 
corresponding expressions to 0(e) at the non trivial fixed point in the O(n) <p theory 
which is found in the e = 4 — d expansion around d = 4. Restricting to the two point 
function of T nn , where n denotes the normal component on the boundary, we are able to 
recover the results of Eisenriegler et al. [4] . 

In section 4 we further consider general results for the energy momentum tensor as 
defined by variations with respect to the metric g^ v on a curved manifold with a smooth 
boundary when variations of the induced metric, and also related geometric quantities 
on the boundary, are taken into account. The implications of diffeomorphism invariance, 
which now go beyond the usual conservation equation V^T^ U = 0, and Weyl invariance 
under local rescalings of the metric, which extend the traceless condition g^ u T^ = 0, 
are derived. These conditions relate components of the energy momentum tensor on the 
boundary to new local boundary operators present in any field theory defined for arbi- 
trary smooth boundaries. The results are verified for a general scalar field theory treated 
classically or to zeroth order in the loop expansion. For completeness a brief summary of 
the essential results used here arising from a covariant geometrical treatment for smoothly 
curved boundary surfaces is given in appendix A. Appendix B contains a discussion of the 
extension of the treatment in section 4 to fermion fields while appendix C contains some 
of the salient details of the 0(e) calculations necessary to obtain the results in section 3. 

2 Conformal Invariance with a Plane Boundary 

For a flat <i-dimensional space with coordinates x^ = (xi, x) we assume a plane bound- 
ary at x\ = 0. It is then necessary to restrict the conformal group to the subgroup leaving 
x\ = invariant [2]. Besides d—1 dimensional translations, 0(d— 1) rotations and scale 
transformations 

Xi * Xi ~t~ (Xi , Xi RijXj , X ^ ^ \x ^ , (2.1) 

where i = 2, 3, . . ., this restriction also allows special conformal transformations 

Xfl X "n(l) X2 ' n (*0 = 1 + 2 *-* + ^a; 2 . ( 2 -2) 
so long as b\ = 0. In this case it is easy to see that for two points x^, x'^ 

( _ '\ 2 (x — X ) 2 x± , x'i , . 

[x x ' ^ n(x)n(x') ' Xl ~" n(x) ' Xl ~* n( x >) ' { ' 

so that it is straightforward to form an invariant under the restricted conformal group 
which we take as ^ 

V = (x-x'^ + Ori + x;) 2 ' ( ' j 
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In consequence for a scalar operator 0(x) of dimension rj the connected two point function 
has the general form, as required by conformal invariance near a boundary, 

(0(x)0(x'))= ^— L_- F(i;), (2.5) 

so that for x\ = or = the magnitude is given by F(l) (if O(x) obeys Dirichlet 
boundary conditions F{1) = 0) whereas the scale of the bulk amplitude neglecting surface 
effects is given by F(0). 

For the energy momentum tensor T^ v under a general conformal transformation x — > x 

doc 

T^ v {x) -> 0(x) d ^ Ma (x)'^. I//3 (x)T Q/ 3(x) , TZ^(x) = Q.(x)-^- , K^n^ = 8^ . (2.6) 

It is convenient to write the two point function for the energy momentum tensor in the 
form (in this semi-infinite geometry (T^ v ) = 0) 

{Tni/(x)T ap (x )) = - Tyzd ^-m^°"p( s ' x i) ' s = x — x . (2-7) 
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figure 1 

By a translation and an appropriate special conformal transformation of the form 
(2.2), with bi = 0, we may choose x = x' = so that x, x' lie on a perpendicular to the 
boundary as in fig. 1. In this case, by invariance under rotations and scale transformations 
as in (2.1) which preserve this perpendicular geometry, we may write 

Aim = a , Aijii = (38ij, A llk i = Ski , Anki =n /Sik, ^ 
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with other components given by symmetry (Ap Uap = A vpap = Ap Vp(J ). The coefficients 
a, P, j3' ,7, 5,e are functions of 

v = , Xl =y>x' 1 =y' . (2.9) 

y + y' 

Assuming the form (2.8) it is trivial to obtain the conditions for tracelessness 

a+(d-l)p = 0, (3 = 0' , (2.10a) 
(3 + (d- 1)5 + 2e = . (2.106) 

In addition it is necessary to impose the conservation equation dp{Tp^T ap ) = 0. To 
derive these extra conditions* it is necessary to extend the form (2.8), valid for s = 0, to 
s ^ 0. For an infinitesimal transformation leaving x'^ = (y', 0) invariant 

5xi = b{X 2 — 2xib-x. — biy' 2 , 5x\ = — 2xib-x , 

%',o)=o, 5( y ,o) = (o,o, £ = Hy 2 - y' 2 ) , 



then the 0(d) transformation matrices as defined by (2.6) at x and x' respectively are 

n ( 1 -Ujv\ n' -( 1 -W\ (oi9) 



Hence 



A ni i(t y, y') = -^-^ £i (ya -y(3- 2y' 7 ) + 0(£ 2 ) (2.13) 

v - y 



and the condition d^TnTn) + <9i (TuTn) = gives, using (2.10a), 

(y 2 -y' 2 )d y (3 = 2y'(d(3- 1 ) . (2.14) 

Similarly 

2 

Aiikiit y, y') = -o a ($ke&(yP - y5) - (6 ik £t + 5 u £ k )(ye - y' 7 )) + 0(£ 2 ) , 

V / (2.15) 

AiMt, y, y) = 2 _ , 2 {Sijtkb/P - y'S) + (5 jk & + 5 ik ^){ yi - y'e)) + 0(e) , 
y y 

and the associated conservation equations, using (2.10a,b), lead to (2.14) again and also 

(y 2 - y' 2 )d yl = 2y'(d 7 - (3 + 5 + de) . (2.16) 



* The following arguments are an adaptation of a similar discussion of conformal invariance require- 
ments on three point functions for general d [5]. 



(2.14) and (2.16) exhaust the restrictions following from restricted conformal invariance 
and the conservation equation for T^ v in the perpendicular geometry assumed for (2.8). 
In terms of the variable v in (2.9) (2.14) and (2.16) simplify to 

v jLp = d/3 - 2 7 , v-^--f = d-f- /3 + 5 + de . (2.17) 
dv dv 

Even after using (2.10a,b) to eliminate a and, for instance, e there remain two coupled 
linear differential equations amongst three unknowns so in general an arbitrary function 
of v appears in the solution. This is exemplified by particular cases subsequently. At large 
distances from the boundary the solutions should tend smoothly to 

a(0) = (d-l)C, (3(0) = 8(0) = -C , - 7 (0) = e(0) = \dC . (2.18) 

For d = 2 C is proportional to the Virasoro central charge while for d = 4 C is related to 
the coefficient of the Weyl tensor squared in the trace of the energy momentum tensor on 
a curved space background [6]. 

However for d = 2 the index i is restricted to just i = 2 and in (2.8) A^ uap depends 
only on the combination 5 + 2e. This is reflected in eqs. (2.10b) and (2.17) so that they 
can be solved to give 

a = _P = S + 2e = C(l + v 4 ) , 7 = -C(l-f 4 ). (2.19) 

The form of 7 corresponds to T i2 = on the boundary x\ = 0. 

As a further application we may consider the correlation function of the energy mo- 
mentum tensor with a scalar field O of dimension r\. In this case as well as (2.6) for 
x — > x 

0(x) -> fl(x) v O(x) 

and instead of (2.7) 

(T^(x)0(x')) = {x _ l x , )d+ri S^(s, Xl , x[) . (2.20) 

In the perpendicular geometry x^ = (y, 0), x'^ = (y' ', 0), y > y', then imposing traceless- 
ness gives 

S n = w(v) , Sij = -j^j w(v) S i:i . (2.21) 
Following a similar discussion to before 

Suit y, y') = 2^ w(v) + O(e) 

y z — y' z d — 1 

and the conservation equation diTu + d{Tn = leads to 

d , N d — T) r) _ 

—w(v) H -w--w = 0. 2.22 

dv 1 — v v 



This has the solution 

w(v) =Sv ri {l-v) d - ri , y>y'. (2.23) 
Conversely if y' > y then, with now v = (y' — y)/(y + y'), the solution becomes 

w (v) = Sv r, {l + v) d - TI . (2.24) 

The coefficient S appearing in (2.23,24) may be related to (G(x)) by using the operator 
product expansion of T^ v {x) and O(x') for x — > x' which with our normalisations has the 
form 

T^(y, 0)O(y', 0) ~ ^ 0(y', 0) , 
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(2.25) 

A 11 = -A u , ^ = --<%, S d = _^. 
Hence from (2.20,21) and (2.23) or (2.24) with v -> 

<°("' » = -T(2?F- (2 ' 26) 

Cardy [3] derived further relations by considering analogous short distance expansions 
in the neighbourhood of the boundary where 0(x) is expanded in terms of low dimension 
boundary operators. Supposing that the expansion is restricted to the unit operator and 
the energy momentum operator T^ v itself on the boundary (where Tn = 0) he assumed 

C%0)~ {O(y,0))(l + by d T 11 (0,0) + ...) , (2.27) 

then from the solution (2.25) for y' — > and also (2.7,8) 

2/ 2d o 

This shows that the coemcient b does not depend on the particular operator O. In two 
dimensions from the solutions (2.19) a(l) = 2C and this leads to so called hyperscaling 
relations [3] for b which is then independent of particular boundary conditions and which 
appears in universal terms in the expression for the Casimir energy for parallel plate 
geometries. For d > 2 there appears to be no general relation between ct(0) and a(l). 

Alternatively if in (2.20) with the solution (2.24) we take y — > then in this case the 
limit is non singular and from (2.26) we obtain 

S'dy 1 



(T 11 (y,O)Tn(0,O)> = a(l) 7 2 5 , a(l) = —2 d — . (2.28) 



(T n (0, 0)O(y', 0)) = -^ (0(y', 0)) . (2.29) 



More general configurations than the perpendicular geometry discussed above may 
be obtained by conformal transformation. If x'^ = (y', 0) is fixed then we may take, 
generalising the infinitesimal transformation (2.11), 

1 + b 2 y' 2 v 2 - v' 2 

g,,o) -(*!,*:), ^ = Y^ y ^ x = rw b (2 - 30) 
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and 

/ 1-&V 2b jy 

1+b 2 y 2 1+b 2 y 2 

V l+b 2 y 2 °V l+b 2 y 2 

and correspondingly for 1Z' with y — > y'. As the magnitude of b varies x M moves on a 
semi-circle from (y, 0) to (y /2 /y, 0) for 6 2 — > oo, as shown in fig. 1, and v — > —v. It is easy 
to verify from (2.30) that (y — y') 2 — > (x — x') 2 and v 2 = (y — y') 2 / (y + y') 2 — > v 2 as given 
by (2.4). Hence (T pu (x)T ap (x')) may be found for arbitrary x,x' in terms of a, (3, 7, 5, e. 
Assuming (2.6) and (2.7,8) with (2.30,31), taking x[ = y', gives 

(T^(x)T ap (x')) = (S(v) 6p V 6 ap + e(v) (l^(s)I up (s) + I pp (s)I U(7 (s)) 

+ (P(v) - S(v)) (X^X,, 5 ap + X' a X' p 5^) 
- ( 7 ( V ) + e(v)) {X^X' a I vp {s) + X v X' p 

+ X^X' p I V(7 (s)+X u X' a I pip (s)) 

+ (a(v) - 2(3(v) + 4ry(v) + 6(v) + 2e(v))X^X v X' a X' p 



where v is given by (2.4), s = x — x' and 

S 



(2.32) 



= n^ a n a = N(x? - x[ 2 - s 2 , 2x 1 s) , 
X' a = - K' aa n a = I^{s)X p = N(x[ 2 - xl - s 2 , -2x[ s) 
iV- 2 = s 2 ( s 2 + 4a: 1 x / 1 ) = ( S 2 ) 2 /z; 2 , 



(2.33) 



for n a = (1,0) defining the normal to the boundary. It is easy to check that (2.32) is 
consistent with T M/i = 0, using (2.10a,b), since X, X' are unit vectors. At large distances 
from the boundary, when v — > 0, all terms containing X or X' vanish by virtue of (2.18) 
and the usual result [6] for no boundary is obtained. 

As a special case we may obtain 

(T 11 (0,x)T 11 (0,x / ))=«(l) { ^_\ f)2d , (2-34) 
which may also be found more directly from (2.28). In a similar fashion from (2.29) 

(T n (0,OK%,x)> = -± (^-) W*)> . (2.35) 

3 Calculations in Specific Models 

The simplest conformal field theory for arbitrary d is perhaps that corresponding to 
a free scalar field cf> for which the energy momentum tensor is 

T„ v = d^d v 4> - \ ((d - 2)d„d v + 5^d 2 )(P 2 . (3.1) 
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This is conserved and traceless on the equations of motion d 2 (p = 0. In order to maintain 
conformal invariance with a plane boundary at x\ = it is sufficient to impose Neumann 
<9i0(O, 0) = or Dirichlet 0(0, 0) = boundary conditions. In the perpendicular geometry 
described in the previous section the basic two point function is 

my,0)M = ^^( ¥ ^± w ^) (3.2) 
and we may easily find 

{Mb, 0)<W, 0)) = - {didjtty, 0)<P(y', 0)> = i- (^rp ± (y + y/jO ' 
(did j( l>{y, ti)d k d t <i>{tf , 0)> = £ (SijSkt + S ik S je + 5 ie 5 jk ) ( ^_* /|d+2 ± ( y + y /)d+ 2 ) • 
With these results it is a matter of straightforward calculation to find 



(3.3) 



Sja(v) = l + v 2d ±l(d-2)d^ v d ~ 2 (l - v 2 ) 2 , 

Sh(v)= -^TTJ^-^i v d ~\l-v*)) , (3.4) 

Sje(v) = (1 + v 2d ) ± \ {(d - 2)(v d ~ 2 + v d + 2 ) + 2dv d ) , 

with (3,8 determined from (2.10a,b). It is easy to verify that these results satisfy the 
differential equations (2.17). For v -> (3.4) gives C = (d - l)~ 1 S d ~ 2 in (2.18). 

For the scalar operator (p 2 , of dimension d — 2 in this free theory, then from (3.2) after 
subtraction of the usual short distance divergence 

{M) = ± J^W d lM^- (3 ' 5) 

and it is easy to verify that this is compatible with (T^ u (y, 0)(j) 2 (y', 0)) as determined from 
(2.20,21) and (2.23) or (2.24). 

An equally simple conformal field theory is that provided by free fermion fields for 
which the energy momentum tensor is 

T^u = §$(7/X + 7^)V> ■ (3-6) 
The basic two point function is taken as 

= ^ C ( ' {X ~fJ + U Zf'^} ) , ^ = (-xi, x) (3-7) 
bd V [x — x ) a {x — x ) a / 

where U is a matrix satisfying 

t/ 7l = - lx U , t/ 7i = 7i t7 , t/ 2 = U 2 = 1 . (3.8) 



This corresponds to boundary conditions 

( 1 -Wl^= ' ^ 1 - lJ )\ dM =0. (3.9) 
It is straightforward to show that in this case 

e(v) = -\d8(v) , p(v) = 8(v) = -^(l+v 2d ) , 7 ( v ) = -Jd^(l-t; M ), (3.10) 

which is similar to the scalar field case after averaging over Neumann and Dirichlet bound- 
ary conditions. 

For a less trivial calculation of the universal functions defining the two point function 
of the energy momentum tensor we consider the well known critical point in scalar field 
theories with a 4> 4 interaction where g* = 0(e). At this non Gaussian fixed point the theory 
is conformally invariant and moreover critical exponents and other universal quantities may 
be calculated as an expansion in e = 4 — d using standard perturbative loop expansions. 
Thus for an n-component field with an 0(n) invariant interaction 

v($) = ±g(4> 2 ) 2 ( 3 - n ) 

we consider the 0(g) contributions to e and 7. These involve only off-diagonal components 
of T^v so that the expression (3.1) is still sufficient. The relevant Feynman diagrams 
are calculated in appendix C for either Neumann or Dirichlet boundary conditions. The 
leading terms in the e expansion are given by (3.4) while the order 0(e) corrections from 
(C.7a,b) and (C.9a,b) are 

e^(v) = Ke \ tv 2 (1 + 4v 2 + v 4 ) In (1 ~f ^ T f v 2 (l + v 4 ) =F 2v 4 



+ (-3(1 + v 8 ) + ±v 2 (l - v 2 ) 2 )ln (1 J Y 

- 3(1 - v 8 ) Inv 2 + \v 2 (l - v 4 ) (l - 0^ + j^y^ 2 
-5v 2 (l + v 4 )-lv 4 + f^^Y (3.12a) 



,2\2 

1 

>:■" 

\\nv 



W (v) = Ke I ± hv 2 (l - v 4 ) (in ^— ^ - i) 

+ ( 3( l_^ ) _ t; 2 (1 _ w 4 ))ln (lz : Z!) 2 



+ (3(1 + v 8 ) -v 2 (l + v 4 ))lnv 2 



where 

In n + 2 

K -Sj9^T8- (3 - 13) 

It is straightforward to verify that these results satisfy (2.17), together with (2.10a,b), 
which in this context become 

(yA. _ 4 ) 7 « = - | /?« + ip e (D , - A)^ = -2 7 « . (3.14) 

From these equations 

pW(v) = Ke | ± 5^ 2 (1 - z; 2 ) 2 In (1 ~f )2 =f fv 2 (l + v 4 ) ± 10t; 4 



+ i(l-f 2 ) 2 (3 + 4t; 2 + 3?; 4 )ln 



2\2 



(3.15) 



Inv' 



+ m-v*)\nv 2 } T 

z 1 — v 2 

-2v 2 (l + v 4 )+2v 4 



All quantities 7^, Z^ 1 ) vanish as i> — > since there is no change in C, given by (2.18), 
at this order. As required by the boundary condition Tn = for x\ = 7 < - 1 - ) (l) = 0. 
However e^(v) is singular as v — > 1 although /3^ 1 - ) (t') is well behaved in this limit. In the 
next section it is shown that Tn = — is a well defined operator on the boundary x\ = 
but this does not restrict the traceless part of to be finite in general on the boundary. 
Of course for d = 2 all components of T^ v are non singular near the boundary. From the 
above expressions for (3 for v = 1 it is easy to see that the leading terms in the e expansion 
give 

n / , 5 n + 2 



which coincides exactly with the results of Eisenriegler et al. [4] who calculated the two 
point function of Tn at x\ = as in (2.34). 

4 Energy Momentum Tensor on Curved Manifolds with Boundary 

In a quantum field theory defined on a manifold Ai with an arbitrary smooth metric 
g^v the energy momentum tensor may be defined as a finite local composite operator by 
considering variations with respect to g^ v . For E M then if M, of dimension d, has 
a boundary dAi, of dimension d — 1, this may be determined by a^b(x) for x arbitrary 
coordinates on dAi. On <9.M the induced metric is defined by 

7y(*) = ^(x b )e^(x)e^(x) , e^(x) = ^ , (4.1) 
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and we may also define related geometric quantities such as the unit inward normal n M (x), 
satisfying n M e M i = 0, n^n^ = 1, and the extrinsic curvature -?Qj(x) given by (A.l) which 
has dimension 1. 

For any quantum field theory the vacuum energy functional W(g,Xh), which here 
depends on the metric g uv and the embedding of dM in M given by x^, may be defined 
by a functional integral of the generic form 

e w = Jd[(j)}e- S ^ (4.2) 

for fields (p on M satisfying suitable boundary conditions on dM and So the bare action 
including all necessary counterterms in an appropriate regularisation scheme so that W 
is a finite functional for arbitrary smooth metrics g nv and boundaries x£. Recently [7] 
we have shown how this may be achieved to two loops using dimensional regularisation in 
four dimensions for renormalisable scalar field theories with quantum fields obeying either 
Dirichlet or generalised Neumann boundary conditions. 

In such a framework we may write under variations in the metric 

-5 g W= [ dv\8g^(T^) 
Jm 

+ [ dS(±6f*{B ij ) + 6n' i {\ li ) + %6K ij {C i *) + ...) , 

JdM 

where dv = d d x^fg, dS = d d ~ 1 x\^y. This relation defines insertions of local operators 
Bij (x), A n (x), Ai(x) and C lJ (x) on dM. The neglected surface terms in (4.3) involve vari- 
ations of higher dimension geometric tensors on dM, such as components of the curvature 
tensor like R n inj = n ll e v in a e p j R^ vap -, but these are absent in simple theories. 

To derive conservation equations we assume invariance under diffeomorphisms, or 
reparameterisations of the coordinates, as given by 5 v x fl = —v fl (x). On the boundary 
dM we also require 5 v x l = — ?/(x) where v^(xb) = e M j(x)?/ (x) + n Ai (x)w n (x). On M the 
external metric is required to transform as 

&v9nv = + V v v u . (4.4) 

The corresponding variations on dM in 7^, n M and induced by 5 v g av in (4.4) are 
obtained in the appendix to be 

Svlij = ViVj + VjVi - 2Kij v n , 
8 v n» = - V n ^ - e^idiVn + K ijV j ) , 

^v-^ij — -\- [Rninj K{ k Kj jl< n -(- V \djV n j 

C v Kij = v k V k K XJ + tiV k K kj + VjV k K ik , 

for Vj the covariant derivative acting on tensor fields over dM with the Christofiel connec- 
tion prescribed by the metric 7^. For a general difieomorphism it is necessary to assume 
also a shift in the boundary surface represented by 8 v x£(x) = — n M (x)v n (x). 
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(4.3) 



(4.5) 



By integration by parts invariance under diffeomorphisms on M. gives the usual con- 
servation equation 

V^(T^) = 0, (4.6) 

(for simplicity we neglect contributions from any sources for other local operators on which 
W may also depend). On the boundary invariance gives from the results (4.5) 

\ . . (4.7) 

(T ni )= -V j (Bij) - \ViK jk {Ci k ) +Vj(K ik (Ci k )) . 

In addition, after setting (A M ) = 0, we may also obtain 

- n *TT W = ( T nn) + K^iBij) + \ (R ninj - K lk K k ) (C*) + JViV^C*') . (4.8) 

(4.7) and (4.8) ensure that T nn = n^n v T^ and T ni = n^e v iT^ u are well defined local 
operators on dM. 

If we also assume Weyl invariance under local rescalings of the metric 

5 a g^ = 2ag^ (4.9) 

then on M this requires as usual 

g^(T^) = 0. (4.10) 
From (4.9) the induced variations on the boundary are 

5^ = 2a Y J , S^n" = an» , 5 a K tJ = a K tJ + %d n a . (4.11) 

Weyl invariance then also requires 



fi(B ij ) + ±K ij (C i i) + (\ n ) = 0, 
%{C*) = 0. 



(4.12) 



For d = 2 and a plane boundary x\ = (4.7) is equivalent to the operator relation on the 
boundary T\ x = —d x B, for x^ = x, while from (4.12) for conformal invariance B = 0*. 

As an illustration of these results we a consider a simple scalar field theory specified 
by bulk and surface contributions to the action of the form 

S(4>) = I dvC , £= Id^d^cj) +\tR4> 2 + V{(f>) , 

J f ' ' (4.13) 

S(<p)= / dS (Q(0) + \ P K <p 2 ) , 

J dM 



* Similar equations have been obtained by Cardy [8] . 
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where R is the scalar curvature on M. and K = ^Kij. In the corresponding quantum 
field theory to lowest order in the loop expansion we may write 

W(g,x h )^ = -S(<P)-S(<j ) ) , (4.14) 

where is a solution of the classical field equations and boundary conditions such that 
S + S is stationary, 

V 2 - T Rcf> - V'{cj>) = , («9 n - pK <j> - Q'(<l>))\ dM = . (4.15) 

Under a variation in the metric, since 5 g R = 5g^ u R fll , — (V M V^ — g^ u V 2 )5g^ , 

S g S(cf>)= [ dv\8g^T^ 
Jm 

+ [ dS ±r((- h^VnSg^ + n^eJVrfgn^ 2 + h^dg^d^ 2 - n^dg^d^ 2 ) , 

JdM 

(4.16) 

for = g^ v - n^n u and 

Tfjw = d^d^ + tR^ cf 2 - r(V M V, - g^V 2 )^ 2 -g^C. (4.17) 

It is easy to check that V M T Mi , = using the equation of motion (4.15). For V = and 
r = \{d — 2)/(d — 1) it is also straightforward to verify that g^ v T^ v = on flat space 
and that T Mi , is then equivalent to that assumed in (3.1) for the free scalar conformal field 
theory. Using (4.16) and the results in the appendix it is possible to write 

8 g (S + S)= [ dv\5g^T^+ [ dS {\8f j B i:j + <^A M + ^SK^C') , (4.18) 
Jm JdM 

where, using the boundary conditions on (j) in (4.15), 

B tj = % (- W) + 2rQ'((f>)(f> + 2rpK<p 2 - \ p K<f) + (p - r)K tj <j> 2 , 
& = (p - 2r)f V , A M = 0. 

The general results derived above may now be verified directly. On the boundary from 
(4.17) 

T ni = d n cf>d i( f> + rR ni <P 2 -t d z d n <p 2 - t K?d j( j> 2 , (4.20) 

where d n (f) may be eliminated from (4.15). Using the Gauss-Codazzi equation R n i = 
diK — V ' jKj it is then easy to show that 

T ni = -V j B id - \ViK jk C> k + Vj(K ik Ci k ) , (4.21) 

which is in accord with (4.7) as expected. On the boundary also 

T nn = d n cf>d n( f> - - \tR$ 2 - V{<f>) + rR nn( p 2 + tV V - TKd n( p 2 , (4.22) 
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and using (A. 7, 8, 9) we may verify (4.8) in this case as well 

n^(S + S)= T nn + K^Bij + \ (R ninj - K lk K 3 k )C^ + \ViVjC^ . (4.23) 

For Weyl invariance from the second of the conditions in (4.12) 7ijC* J = we obtain 
p = 2t and then 

y'Bij = (d - l)(-Q((P) + 2tQ'((P)<P) + r(4(d - l)r - (d - 2))K(P 2 . (4.24) 

Hence 7 lJ Bij = gives the same result for r as g^T^ = 0. If Q(4>) = \ccj) 2 the remaining 
condition is satisfied if either c = 0, and d n (t>\dM = 0, or if c — > oo, c</> = 0(1), so that 
(f>\dM = 0, showing how both Neumann and Dirichlet boundary conditions on scalar fields 
are separately compatible with conformal invariance. 

5 Conclusion 

It is clear from the results of this paper that two-point correlation functions of the 
energy momentum tensor in conformal field theories with appropriate boundary conditions 
preserving conformal invariance may have a dependence on the conformal invariant variable 
v, defined here by (2.4), which depends on the particular conformally invariant theory and 
its boundary conditions whenever d > 2. An interesting question, beyond the scope of our 
considerations here, is whether such functions are measurable in realistic statistical physics 
models at a critical point when we may take d = 3. For models in the same universality 
class as the Ising model then our s expansion results would perhaps be relevant setting 
e = 1. For such applications it would be desirable to find functions which agreed with our 
results to O(e) but remained solutions of the equations (2.10a,b) and (2.17) for general d 
since they would then extrapolate to the unique functional form for d = 2 given by (2.19). 

We are grateful to John Cardy for stimulating conversations and sending us a copy of 
ref. 4. 
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Appendix A 



For application in section 4 we here summarise the essential results in a geometrical 
treatment of a boundary dM, parameter ised by coordinates x l , where dM is specified in 
terms of the coordinates for M by a^(x). A natural tangent frame basis on dM is given by 
e M i(x), n^(x), with n M the unit inward normal, and the extrinsic metric on dM is defined 
as in (4.1). The symmetric tensor Ki 3 forming the extrinsic curvature and a connection 
Y k 3 are defined on dM by 

V, e M. = Q^. + T^eVi - e%Y% = n»K l3 , 

= d^ + T» p eV = -K ije M , d t = , 

where T^ p the usual Christoffel connection formed from g^ v . Acting on tensors on dM 
a covariant derivative Vj may be defined by the connection = T k t and from (A.l) 

and (4.1) Vi'jjk = so that T k 3 is just the Christoffel connection formed from 7^. From 
(A.l) it is straightforward to derive the Gauss-Codazzi equations relating the Riemann 
curvature tensor R^ vap for x G dM, with zero and one component along the normal n M , 
to the intrinsic Riemann curvature Rijki of dM associated with the covariant derivative 
Vj and also the extrinsic curvature K i3 . 

To derive the implications for 7^, n M and of the variation of the basic metric (4.5) 



induced by a diffeomorphism we first note that from (A.l) e^e^V^ = Vtv 3 — K 
which leads directly from the definition (4.1) to 5 v ^i 3 in (4.6). The variation of the normal 
vector in (4.6) may be determined from the equations 5 v n p e^i = and 5„n M n /i = 
— \S v g li ' v n p n v = n^n^V \ 1 v v . For the variation in the extrinsic curvature Ki 3 = n^Wie^j 
from (A.l) in order to obtain the result in (4.6) we use 

8 v Kij = n^n u V^v v K i3 + n p 5 v T£ e a i ef > i , 



crp 

5 v T£ p = g^V w V p) v u - v x R x[ap f , R k{ij)n = V {l K j)k - V k K, 



IJ 7 



n^ePjVaVpv" = - K l3 n^n v V p v u + V l d j v n - K lk K k v n + V 3 v k K ik + V l (K jk v k ) . 

To express the variations in the metric and its derivatives on dM in terms of variations 
in 7y, n M and K i3 , as required for the simplification of (4.16), we use 

Sg^ = 5n^ n v + n» 8n" + e^ 3 6f j . (A.2) 

It is easy to see that, for h pv = g^ — n p n v , 

KJg^ = %5fi , n^eJSg^ = eJSn" , (A.3) 

and also, using (A.l), 

n^eJViSg^ = Vi{ej8n v ) - 2Kn p 5n» + K i3 5f j . (AA) 



15 



To determine a corresponding expression for h^V n Sg^ v we first use K i3 - = —e^ie^jV^n^ 
to find 

SKij = - 5n v + e^ 3 n x 6T* V 

= - e^e^V^ Sn v - e^V^Sn^ + K tj n f ,Sn t " - % (l K j)e S^ ke + \e^e vj V n 8g^ v . 

From this it is straightforward to obtain 

h^VJg^ = 2SK - 2K n M 5n M + 2V { {ejSn*) . ( A5) 

Finally we consider the variations induced by a shift in the boundary surface along a 
normal as given by S t x£ = —n^St. We let 



8' t ePi = diS t x^ - 6tn a T£ v e"i = StK i3 e™ - n^d.St , 
S' t n^ = dtn^ - n v Y v afJL n a St = eJdiSt . 



(A6) 



With the induced metric given by (4.1) we have 

S t % = 2K l3 St , (A7) 

and using = \e tl j from (A.l) 



S t Kij = e k d k St Vie" j - n M [V n , V,]e^ + n.V^Je^ 
= - VidjSt - (R nj m - K jk K k i)5t . 



(A8) 

For integrals over a local scalar function / on M or, restricted to the boundary, on dM 



S t dvf = dS Stf\dM , 
'M JdM 



St I dSf\ dM = [ dS(St(-d n f + Kf) + S t f) 

JdM JdM 



(A.9) 



dM 



Appendix B 

The discussion in section 4 of the energy momentum tensor and related boundary 
operators, as defined through variations of the metric for a curved space background, 
seemingly requires modification if fermion fields are present. As is well known [9] it is 
then necessary for a consistent treatment to introduce a tangent frame basis V^ a such 
that g^ v = V^aV^a and a corresponding connection u^ a b = —^^ba defined by V lx V v a = 
dpV v a + V v iia V a a + uj^abV^b = 0. The conventional Dirac action is 

S D = [ dvi;{^% + M)V , r = V^ala , = c> M + \uj^ ab a ab , (B.l) 

J M 
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for "f a the usual Dirac matrices {7 a , 7&} — 25 a ^, a a b — ^[7a> 

7 6 ]. On the classical field 
equations, (7 M V^ + M)ip = 0, ^(7^^ — M) = 0, then S D = 0. 

When considering variations in F M a , which gives a variation in the metric g^ v ', it is 
also necessary to allow for contributions proportional to 

^ = 1 {5 y» a y» a _ 6V » a y» a) . ( S .2) 

By straightforward calculation, using Su^ab = K-[aV /i 5V I &] — V^oV^Vt,^^, assuming 
-0, ^ satisfy the equations of motion but without any specific boundary conditions then 

S V S D = [ dv(\5g» v T°+\^i}[a^M]i))-t dS\t^ ln a^ , 

Jm ' ~ (5.3) 

As usual [9], assuming invariance under local rotations which in this context requires 
[<Jab, M] = 0, terms proportional to in SvS D on M are absent. Clearly from (B.3), 
for the simple Dirac action S D , there are no contributions under a variation of the metric 
proportional to ... as in (4.3) or (4.18). 

Under variations in ip, ip, subject to the field equations on M, 

S^S D =[ dS{8$ ln il>-$ ln 8il>) . (5.4) 

JdM 

Conventionally [10] linear boundary conditions on ip',ip are chosen, as in (3.9), so that 

Vln^\ dM = , (5-5) 
and hence i7 M V^ is a symmetric operator and in (B.4) we also have 8^S D = 0. 

The basic equations derived in section 4 from invariance under diffeomorphisms giving 
expressions for T n i , T nn on dM. in terms of boundary operators now require extension to 
take account of terms proportional to t^ v as well as terms involving SY J \ ■ ■ • from variations 
of the metric. Corresponding to a diffeomorphism, leading to (4.4), we require 

5 v V» a = v u d a V» a - d a v» V\ = -V ff u" V\ - v a uj aah V% , 

and hence from the definition (B.2) 

V? = 1 ( W - W) + v a u aah V\V v h . (B.7) 

The essential identity of section 4 has the form in this case 

S V S D = - [ dS (v n T° n + v l T°) 

JdM 

+ § / dS (v^abm^aab}^ + VV'^7 n cx^) (5.8) 

JdM 

+ \ \ dS (6 v i}ry n i() - V>7rAV0 • 
JdM 
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(5.6) 



This may be verified using Syip^dyip from (B.6) and from (B.3) 

(B.9) 



Since S D = there are no terms involving n^SS /5x£ as in (4.23). To derive (B.9) it 
is useful to note that we may write the tangential components of the spinor covariant 
derivatives on dM. in the form 

Vi = V; + \Ki^ n ^ , Vi = di + u>i , 

(-D.1U) 

Vi7 n = da n + [u>i, 7 n ] = , Vi7j = . 



If the boundary conditions on -0, imply also V'TnO'ijV' = 0; besides (B.5), then the 
additional surface terms in (B.3), (B.8) are zero and = = on dM and the 
discussion in section 4 does not require any modification. 

Appendix C 

Here we outline the essential steps in the calculation of the two point function of the 
energy momentum tensor at the critical point to 0(e). We calculate only those contri- 
butions to e and 7, as defined by (2.8), which involve off diagonal elements of since 
these are not affected by interaction terms in this theory. By virtue of (2.10a,b) and (2.17) 
knowledge of e and 7 is clearly sufficient to determine all other pieces of the the two point 
function for T^ v as well as providing a convenient consistency check. 




figure 2a figure 2b 



The two basic graphs whose contributions need to be calculated to 0(g) are shown 
in fig. 2a,b. Corresponding to fig. 2a it is sufficient to find the one loop corrections to 
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{(j)(x)(p(x')). It is easy to see that 

/ n^\0) - 1 ! 



(<P(y, 0)0(2/', 0)>« = - -(n + 2)g s s {d _ 2) s Hy, y') ± G(y, -y')) , 



— oo 



X 2 =x 2 + ( 2/ -^) 2 , X' 2 = x 2 + (y' - zf . 

Since the critical coupling g* = 0(e) it is sufficient to evaluate G(y,y') for d = 4. The 
integral in (C.l) then has a linear divergence which may be regularised by imposing a 
cut-off \z\ > e and we find 

G(y, v') = - 2 < i* + 1 > . (c . 2) 

0/ + y') {y-y'r e M + |y I 

For the Dirichlet case the divergence as e — > cancels while in the Neumann case it may 
be removed by a surface counterterm in the action oc (p 2 . At the critical point to this order 
g*/l6n 2 = 3e/(n + 8) and the result for (4>(y, Q)(p(y', 0)) is compatible with (2.5), since 
i] = |(<i — 2) + 0(e 2 ), where the function of the conformal invariant v is given by 

S d (d-2)F(v)^ = -l^e(v 2 ln^±Hl-v 2 )) . (C.3) 

This result is equivalent to that obtained by Gompper and Wagner [11]. For use in calcu- 
lating the contributions corresponding to fig. 2a to (Tjj(y, 0)T k i(y', 0)) we also require, as 
in (3.3) for free fields, 



(d i <j>(y,0)d j <j>(y',0)) = - {didjfty, O)0(y',O)) = — -1—^F^v) 

\y y i 

(didj(j){y, 0)d k d e (f)(y' , 0)) = (d lj d k £ + S ik Sj£ + SuS jk ) - -77^4 F 2 (v) , 

I y y 1 



(C.4) 



(C.5) 



where F l7 F 2 are determined in terms of F by 

= 2r/F(t;) - u(l - f 2 )F» , 
F 2 (u) = 277(277 + 2)F(v) - u(l - v 2 )(4r] + 1 + 3v 2 )F'(v) + v 2 (l - v 2 )F"(v) . 

Corresponding to (C.3) we have, taking 77 — > 1, 

S d (d - 2)F 1 (v)< 1 > = - £ ft; 4 In ± ln(l - z, 2 ) + v 2 ± z; 2 ) , 

n + 8 V v z / 

S d (d - 2)F 2 (v)^ = - 4 £±| e ( v 6 In 1^! ± ln(l - z, 2 ) + v 4 + \v 2 ±v 2 ± \v 4 ) . 

(C.6) 
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With these results we may find the contributions arising from fig. 2a to e and 7 just as 
for free fields in section 3. In the case of 7 the required correlation functions involving 
d\4> may be found from (<f)(y, 0)(f)(y', 0)) and (di(f)(y, 0)dj(f)(y', 0)) by differentiation with 
respect to y, y'. Hence we obtain 

ShP = ^ £ ^ Al + ^ + v * ) In llzp! =F ^ 2 (1 + „ 4 ) =F 2, 4 

- 6f 8 In - 6 ln(l - v 2 ) - hv 2 {\ + v 4 ) - v 4 \ , (C.7a) 



^ 9 ri >, n n + 2 i 9, ^ / (1 — f 2 ) 2 1\ 



-6t. 8 ln^-^- + 61n(l-t) 2 )+t. 2 (l-t) 4 )'. . (C.76) 

1 



For the graph in fig. 2b the corresponding contributions relevant for determining e 
and 7 are 

(Ttjiy, 0)T ki (y', O))^ = - + fofy) Jn(rc + 2)«/ (;^) 2 ---^7 fl 

r°° r 1 1 

.2„./2 / j_ „4 / jd-l„ „4 1 1 



—00 



xv V / dzz d xx 



+ terms proportional to dijdke , (C.8a) 
(T 41 (y,0)T fcl (y / ,0)>J 1) = - <5 ifc §n(n + 2)</ ' ^ ' 



1 1 



x yy' / dzz 4 d d 1 xx : 
J — 00 «/ 



x(x 2 + z 2 -y 2 )(a; 2 + z 2 -y /2 ) , (C.86) 



with X, X' as given in (C.l) and X 2 = x 2 + (y + z) 2 , X' 2 = x 2 + (y' + z) 2 . It is clear from 
(C.8) that in the absence of a boundary the one loop graphs for (T flv (x)T ap (x')) are zero 
[12] so there is no change in the bulk coefficient C in (2.18) to this order and hence the 
results of all calculations for e, 7 at this order should vanish as v — > 0. The integrals are 
more tedious to evaluate in this case but there are no singularities requiring regularisation 
even when d = 4. In this case we obtain 

,2AD_nn+2 ll , 2 2 (1 - v 2 ) 2 4 12^ 



W = 7^ ^ ^ <^(1 - v 2 Y In ^— 1-L- - + 



9 n + 8 5 I v ' v 2 (1-v 2 ) 2 
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- f 2 (1 + f 4 ) (l - (1 _^ 2)2 ) In f 2 - f 4 In v 2 j . (C.96) 

We should note that (C.7) and (C.9) are in accord with the symmetry e(v) = f 2d e(w _1 ), 
7(f) = —v 2d 'j(v~ 1 ) with d = 4. It is easy to see that the apparent singularities as v — > 1 
cancel and in fact 7^ = for ?; = 1. 
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